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ABSTRACT 


A  simplified  analytical  theory  of  the  absorption  of 


energetic  magneto spheric  particles  by  an  inert  satellite  is 


developed  for  the  case  in  which  the  radius  b  of  the  satellite 


is  much  less  than  the  equatorial  gyroradius  r  of  the  particle 
■which,  in  turn,  is  very  much  less  than  the  radius  r  or  the 


satellite's  orbit.  This  theory  and  the  observed  particle- 


absorption  signatures  of  Saturn’s  Ring  G  are  used  to  establish 


an  upper  limit  on  shepherding  satellites  associated  with  the 


ring.  The  resulting  upper  limit,  ignoring  the  absorption  of 


the  optically  observed  particulate  matter,  is  b  =  3-3  km 


(accurate  to  a  factor  of  two)  for  a  single  satellite  or 


2 

Z  b^  =  11.2  km  (accurate  to  a  factor  of  four)  for  an  assemblage 
i=i 


of  n  satellites  of  various  radii  b . .  No  shepherding  satellites 

v 


at  Ring  G  have  been  detected  optically. 
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INTRODUCTION 


During  Pioneer  11 's  passage  through  the  inner  radiation 
belts  of  Saturn  on  1  September  1979,  a  number  of  distinctive 
absorption  features  were  observed  in  curves  of  the  radial 
dependence  of  the  intensity  of  energetic,  trapped  particles 
[Van  Allen  et  al.,  1980a,  1980b;  Simpson  et  al.,  1980a,  1980b; 
Fillius  et  al. ,  1980;  Fillius  and  Mcllwain.  1980;  Trainor  et  al., 
1980;  and  McDonald  et  al. ,  1980].  One  such  feature  was  a  well 
defined  inflection  (not  a  dip)  in  curves  of  the  intensity  of 
protons  (Ep  >  80  MeV)  vs.  radial  distance.  In  the  vicinity  of 
this  feature  the  inbound  and  outbound  curves  (at  a  time  lapse 
of  nearly  four  hours  and  at  quite  different  longitudes)  were 
virtually  identical,  suggesting  that  the  feature  was  a  time¬ 
stationary  one.  Van  Allen  et  al.  [1980b]  designated  this 
absorption  feature  1979  S3  and  labeled  it  as  a  "suspected 
[satellite]  but  interpretation  of  [absorption]  signature 
ambiguous".  The  relevant  ambiguity  stems  from  the  failure  of 
such  a  limited  set  of  data  to  distinguish  the  absorption 
macro -signature  of  one  or  more  small  satellites  from  that  of 
a  continuous  circular  ring  of  particulate  matter  [Van  Allen, 


k 


1982].  The  ambiguity  was  resolved  to  some  extent  in  November 


1980  by  Voyager  1  images  [Smith  et  al. ,  1981],  which  revealed  a 
thin  continuous  ring  of  radius  about  170,000  km,  normal  optical 
opacity  «  3  X  10-^,  and  radial  width  about  500  km  (after 
estimated  correction  for  imaging  "smear").  The  Voyager  investi¬ 
gators  designated  it  Ring  G.  A  detailed  interpretation  of  the 
energetic  particle  data  associated  with  Ring  G  has  been  reported 
[Van  Allen,  1983].  The  mean  radius  of  Ring  G  was  found  to  be 
170,200±  900  km,  in  agreement  with  the  optical  value.  Other  con¬ 
clusions  of  this  work,  which  ignored  the  possibility  of  ring- 
associated  satellites,  were  that  "the  particulates  in  Ring  G  have 
an  effective  radius  R  0.035  cm,  an  areal  mass  density 
a  >  1.4  x  10  g  cm"  ,  and  an  areal  number  density 
n  <  8  X  10“^  cm-2." 

In  related  work,  Goldreich  and  Tremaine  [1979,  1982] 
suggested  that  small,  as  yet  discovered,  satellites  might  be 
closely  associated  with  each  of  the  narrow  rings  of  Uranus  and, 
if  so,  would  act  as  agents  (shepherds)  for  dynamically  maintain¬ 
ing  the  observed  distribution  of  particulate  matter.  This 
theoretical  suggestion  received  resounding  confirmation  with  the 
Voyager  1  discovery  [Smith  et  al.,  1981]  of  two  small  satellites 
(1980  S26  and  I960  S27)  that  apparently  act  as  shepherds  for 
Saturn's  Ring  F  [Gehrels  et  al. ,  1980].  More  recently.  Smith 
et  al.  [1986]  have  observed  two  small  satellites  associated 


5 

with  Uranus'  Ring  €.  However,  there  are  many,  many  narrow  rings 
at  Jupiter,  Saturn,  and  Uranus  that  have  no  thus-far  detected 
shepherds  associated  with  them.  Saturn's  Ring  G  is  one  such 
example;  but,  in  this  case,  no  comprehensive  search  in  longitude 
for  small  (radii  10  km),  nearby  satellites  was  made  by  either 
Voyager  1  or  Voyager  2  [B.  A.  Smith,  private  communication]. 
Hence,  there  is  no  definitive  optical  evidence  either  for  or 
against  the  presence  of  such  satellites  at  Ring  G. 

On  the  basis  of  a  crude  comparative  analysis  of  the 
energetic  particle  intensity  at  Ring  G  and  in  the  inner  slot 
reg-on  (caused  by  Janus  and  Epimetheus),  Van  Allen  [1983] 
reported  that  "there  are  no  satellites  having  radii  of  the 
order  of  a  kilometer  or  larger  in  or  near  Ring  G- "  This 
comparison  ignored  all  of  the  optically  observed  particulate 
matter  in  Ring  G  and  attributed  its  entire  absorption  signa¬ 
ture  to  one  or  more  small  satellites,  thus  attempting  to 
establish  extreme  upper  limits. 

The  present  paper  has  the  same  purpose  but  employs  an 
absolute,  rather  than  a  comparative,  analysis  and  probably 


provides  superior  results. 
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rotational  axis  [Smith  et  al.,  1980;  Acuna  and  Ness,  1980;  Acuna 
et  al.,  1980;  and  Connerney  et  al. ,  1984];  that  the  satellite  is 
moving  in  a  prograde  circular  Keplerian  orbit  of  radius  r  in  the 
planet  s  equatorial  plane;  and  that  the  guiding  center  of  the 
particle's  helical  motion  drifts  in  longitude  in  a  circle  of  approx¬ 
imately  thr  same  radius  and  oscillates  in  latitude  between  mirror 
latitudes  ±  With  all  of  these  simplifications,  the  absorption 

problem  can  be  treated  by  a  simple,  analytical  method  having 
reasonable  accuracy.  Remarks  on  the  effects  of  departures 
from  the  simplifying  conditions  will  be  made  later.  Other 
basic  parameters  are  the  kinetic  energy  R  and  species  of  the 
particle,  its  velocity  vector  V,  its  equatorial  pitch  angle 
%  =  arc  cos  (Vq  .  ^/V^)  (where  Bq  is  the  local  magnetic  field 
vector  and  the  zero  subscripts  denote  equatorial  values),  the 
equatorial  values  of  its  gyroradius  r  and  gyroperiod  T  ,  and 

8  g 

its  latitudinal  bounce  period  Tg-  The  relationship  connecting 
2  6  p  ^ 

ao  and  is  sin  =  cos  \m(4-3  cos  \^)  p .  Extensive  use 

is  made  of  the  formulae  of  Thomsen  and  Van  Allen  [I98OJ.  For 

numerical  examples  of  the  Ring  G  problem,  we  assume  that  a 

representative  particle  is  a  proton  E  =  200  MeV  and  that 

P 

a  representative  value  of  a  is  70°.  The  other  relevant 

o 

parameters  then  have  values  as  follows:  r  =  170,200  km; 

Tg  =  9-084  x  10"2  s;  T0  3.085  s ;  rg  =  2,305  km; 
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\  =  9-6°  =  0.l68  rad;  and  V  =  169,7^0  km  s  .  The  angular 

velocity  m  of  the  satellite  relative  to  the  longitude  of  the 

-2  -1 

particle's  guiding  center  is  2.120  x  10  rad  s  and  thus 
uir  =  3,608  Inn  s-^;  and  the  interval  of  time  between 
encounters  of  the  satellite  with  the  gyrating  particle  is 
2rr/tu  =  296  s.  The  overall  results  are  relatively  insensitive 
to  the  value  of  within  the  energy  range  80-300  MeV.  In  the 
course  of  the  analysis,  numerical  values  of  calculated  quantities 
are  given  in  square  brackets  as  an  aid  to  quantitative  under¬ 
standing. 

During  a  particular  encounter,  there  are  two  opportuni¬ 
ties  for  the  particle  to  hit  the  satellite.  The  probability  of 
a  hit  during  each  opportunity  is  now  calculated. 

The  origin  of  the  adopted  Cartesian  coordinate  system  is 
the  point  at  which  the  guiding  center  of  a  particular  gyrating 
particle  pierces  the  equatorial  plane;  the  origin  noves  in 
longitude  at  the  same  average  rate  as  does  the  guiding  center. 

The  X-axis  is  tangent  to  the  satellite's  orbit,  the  Y-axis  is 
directed  to  the  center  of  the  planet,  and  the  Z-axis  is  parallel 
to  the  rotational  axis  of  the  planet.  In  this  frame  of  reference  the 
particle  moves  along  a  helical  path  on  the  surface  of  a  truncated, 
circular  tube  of  magnetic  flux,  whose  axis  at  the  equator  is  along  the 
Z-axis,  whose  length  is  «  2 r and  whose  equatorial  radius  is  r  . 

The  satellite  drifts  parallel  to  th'  X-axis  through  1  hi s  surface 

with  velocity  v  =  our,  v  =0,  and  v  =0  and  with  its  center 
x  y  z 

a 1  y  -  y_  and  z  =0.  A  hit  can  oe-'wr  only  if  |y  I  £  ( r  *  b). 

o  g 


US 


The  X-Y  plane  section  through  the  surface  is  shown  in 
Figures  1(a)  and  1(b).  For  the  case  shown  in  Figure  1(a),  two 
separate  opportunities  for  a  hit  occur  within  the  two  intervals 
of  time  during  which  the  center  of  the  satellite  drifts  from  A 
to  B  (first  to  second  contact)  and  from  C  to  D  (third  to  fourth 
contact).  Let  l  =  AB  =  CD.  Then  from  the  geometry  of  the  figure. 


L 


•  {y<vb)2 


for  o  £  y  £  (r  -b) 
o  g 


(1) 


For  the  case  shown  in  Figure  1(b),  the  two  separate  segments 
AB  and  CD  merge  into  a  single  segment  AD;  denoting  SL  =  AD/2  for 
this  case, 


A  sample,  composite  plot  of  £/r  vs  y  /r  is  shown  in  Figure  2. 

S  o  g 

It  is  noted  that: 


max 


2  Vbr~  at  |yo|  =  (rg  -  b) 


(3) 
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and  that  the  mean  value  of  1, 


rrb 


or^(yo)  dy  _  ^  _ 

<A>  -  -757—  =  ~ 


(M 


with  integration  over  the  range  0  sjy^js  (r^  +  b)  [cf.  Schardt  and 
McDonald,  1983]-  The  center  of  the  satellite  spends  time  At  in 
tranversing  SL  where 


At 


tar 


(5) 


and  the  radius  of  the  circular  patch  of  intersection  between  the 
cylindrical  surface  and  the  satellite  during  At  varies  from  zero 
to  b  and  back  to  zero  with  a  mean  value: 


<b>  =  rrb/U, 


(6) 


11  uu  1 


except  for  yQ  nearly  equa! 

In  the  vicinity  of  the  equator,  the  particle's  motion  may 
be  visualized  as  composed  of  two  components:  a  Z-component 


[58,000  km  s'1  for  aQ  =  70° ], 


v  »  v  cos  a 
z  o 


(7) 
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and  a  circumferential  component: 


V 


<P 


2tt  r 


T 

g 


[159>^00  km  s-'*"] 


(8) 


The  random  probability  that  the  particle  lies  -within  the 
azimuth  required  for  a  hit  during  At  is 


At  l 

P1  T  air  T 

g  g 


(9) 


If  p1  as  calculated  by  (9)  is  less  than  unity  for  all  values  of  y  , 
then 


<P1> 


<  JL  > 
our  T 


g 


TT  b 

icrT 


g 


[9.6  x  lO-5  b, 
with  b  in  km] 


Sy  (9)  for  Amax  =  2  Vb  rg, 


p^  exceeds  unity  if 


(10) 


b  >  12  km. 


However,  even  if  b  exceeds  this  value  by  as  much  as  a  factor  of  2, 
the  mean  value  of  <p^>  as  given  by  (10)  is  only  5  percent  greater 
than  its  properly  calculated  value. 

The  independent  (assuming  no  correlation  between  gyrophase 
and  bounce  phase)  probability  that  the  particle  lies  within 
±  rrb/U  of  the  equator  at  an  arbitrary  moment  is 


'V  cos  a  Tn 
o  o 


<p2> 


V  cos  a 
B  o 


[For  a  =  70°  and  previously  quoted  values  of  V  and  IL, 

O  D 

<p^>  =  1-75  x  10  5  b,  with  b  in  km.] 

A  typical  mean  value  of  At  (for  b  =  5  tan)  is  0.0044  s;  during 
this  time  interval  the  equatorial  projection  of  the  particle's  trace 
in  space  (-when  near  the  equator)  travels  through  an  azimuthal  arc  of 
0.30  radian  or  0.048  of  a  circle.  A  typical  time  lapse  between 
centers  of  the  two  opportunities  is  1.1  s  or  about  12  gyro- 
periods  of  the  particle,  and  therefore  the  probability  for  a 
hit  in  either  opportunity,  insofar  as  it  is  small,  is  independent 
of  that  in  the  other  (except  for  singular  cases).  When  yQ  «  r^ 
the  probabilities  of  a  hit  during  the  two  opportunities  are 
not  statistically  independent  but  for  b  «  r^  this  qualification 
has  a  trivial  effect  on  the  overall  analysis  and  I  take  the 
probability  for  a  hit  on  a  given  encounter  to  be  twice  that  for 
one  of  the  two  opportunities  therein. 


Thus  the  mean  probability  of  a  hit  (averaged  over  y  ) 
during  either  of  the  two  opportunities  is: 


<p>  = 


<PL>  .  <p2> 


<p>  = 


2  ,  2 
TT  b 


uirVT  T_  cos  a 
g  B  o 


For  a  complete  encounter  (two  opportunities,  assumed  independent) 


<  P  >  = 


2tt2  b2 


our  V  T  T_  cos  a 
g  B  o 


Putting  my  representative  numerical  valufes  in  (13) 


<  P  > 


=  3-36  x  IQ-?  b2 


with  b  in  km.  Equations  (13)  and  (14)  are  appropriate  for 


b  <  20  km  and  are  the  principal  result  of  the  foregoing  analysis 


<P>  in  equations  (13)  and  (lU)  may  also  be  interpreted  as 


the  probability  of  a  hit  within  T ,  the  time  interval  between 

E 


encounters.  If  there  ar«  n  satellites  of  radii  b^,  b^,  •  •  •  bn  moving 


in  random  phases  in  the  same  orbit,  then  the  probabilities 


of  a  hit  are  additive  ani 


<  P> 


_  i  i  _ 

3.36  x  10  I  bf 


v  „>  .« .*V.V  ‘*V 


.V’lVV 
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Equation  (15)  shows  that  <P>  per  T  is  proportional  to  the  sum 

L 

of  the  cross-sectional  areas  of  the  n  satellites,  provided  that 
all  are  large  enough  to  absorb  the  specified  particle  in  a  single 
hit  and  small  enough  (b  <  20  km)  to  meet  the  criterion  for 
applicability  of  (13). 

The  number  of  encounters  per  day  is  86,400/T  =  292  in  the 

case  at  hand.  Hence  the  mean  lifetime  t  of  a  representative 
particle  against  absorption  by  a  satellite  of  radius  b  is  given 
by 


T  = 


1.02  X  10 


28 

days  =  —  years 

b 


(16) 


or  for  n  satellites  having  radii  b, ,  b_,  . . .  b  , 

1  2’  n 


T 


28 


n 

I 

i=l 


years. 


In  (14),  (15)/  (l6),  and  (17),  b  and  b^  are  in  km. 


(17) 
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COMMENTS  ON  VALIDITY  AND  ACCURACY 
OF  EQUATION  (13) 

The  derivation  of  equation  (13),  though  novel  in  detail,  is 
akin  to  those  of  previous  authors  [in  particular,  the  derivation  of 
Schardt  and  McDonald,  1963]  but  has,  I  think,  greater  clarity. 
Nonetheless,  it  is  not  exact  even  for  b  ^  10  km.  A  discussion  of 
its  shortcomings  follows. 

By  equation  (ll),  the  value  of  <p  >  becomes  equal  to  1.0 
if 


V  T„  cos  a 
B  o 

For  b  =  10  km,  V  =  169,7^0  km  s  ^  and  T„  =  3-085  s,  this  condition 

B 

corresponds  to  cosaQ  s  6  x  lO-'*  or  aQ  differing  from  90°  by  only 

0.00U°.  Within  this  range  of  aQ,  <  p^  >  must  be  set  equal  to  1.0 

and  <P>  is  very  much  greater  than  it  is  for  values  of  a  that 

o 

differ  substantially  from  90°.  It  is  therefore  clear  that  particles 

having  aQ  near  90°  are  preferentially  absorbed.  The  dependence 

of  <  P (a  ) >  on  a  is  illustrated  by  Table  1  (all  for 
o  o 

E  =  200  MeV) . 

P 


For  an  angular  distribution  of  unidirectional  particle 
intensity  of  the  form  sin11  ct^,  the  weighted  value  of  <  F  >  over 
the  range  0°  £  £  85°  (i.e.,  ignoring  the  short-lived  range 

85°  £  aQ  s  90°)  is  equal  to  <  P  (aQ)>  with  aQ  =  7*+°  to  17°  for 
n  =  1,  2,  and  3;  and  over  the  range  0°  £  £  80°,  the 

weighted  value  of  <P>  is  equal  to  <  P  (a  )>  with  a.^  =  6b°  to 
70°  for  n  =  1,  2,  and  3.  These  calculations  provide  the 
justification  for  adopting  a =  70°  as  a  representative  value 
for  simplified  considerations  which  do,  of  course,  ignore 
refined  features  of  the  absorption  process  (i.e.,  the  angular 
dependence  of  <P  (aQ)>)*  It  is  noted  that  the  absorption 
signatures  of  Ring  G  were  measured  with  omnidirectional 
detectors . 

If  the  magnetic  axis  of  the  planet  is  tilted  wich  respect 
to  its  rotational  axis,  and/or  if  the  orbit  of  the  hypothetical 
satellite  is  inclined  to  the  planet's  equatorial  plane,  there  are 
additional  considerations.  However,  if  the  tilt  and  inclination 
angles  are  small  (<  1°),  then  the  consequent  effects  are  signifi¬ 
cant  only  for  aQ  near  90°  and  are  otherwise  swallowed  up  in  the 
distribution  of  particle  intensity  with  a .  For  Saturn,  the  tilt  of 
the  magnetic  axis  is  1°  [Connerney  et  al.  ,  198^]  and  the  inclina¬ 
tion  of  the  orbits  of  known  inner  satellites  is  also  ^  1°  [Synnott 
et  al.,  1981,  19831-  Further,  if  the  magnetic  dipole  moment  of  the 
planet  is  offset  from  the  geometrical  center  of  the  planet  and/or 
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if  the  hypothetical  satellite  is  moving  in  an  eccentric  orbit, 


ft 


•Y> 

A 


li 


the  effect  is  to  reduce  the  probability  of  a  hit  ner  T  for  a 

E 

specified  particle.  The  ratio  r  Jr  («  0.0135)  gives  a  measure 

of  the  magnitude  of  the  combined  eccentricities  at  -which  "his 

effect  becomes  significant.  By  the  references  just  cited,  the 

offset  of  the  dipole  is  <  600  km  (i.e.,  <  0.26  r  )  and  the 

8 

eccentricities  of  the  orbits  of  the  known  inner  satellites  are 

^  0.009-  Hence,  the  potential  defects  in  the  simplified  derivation 

of  equation  (13)  as  listed  in  this  paragraph  are  probably  trivial. 

Observational  support  for  this  conclusion  is  provided  by  Figures  U 

and  5  of  Van  Allen  [1983],  which  show  that  the  radial  width  of  the 

absorption  signature  of  Ring  G  is,  in  fact,  not  significantly 

greater  than  2  r  .  The  guiding  centers  of  particles  do,  of 
8 

course,  diffuse  back  and  forth  in  r;  if  the  characteristic  time 
constant  for  such  diffusion  is  less  than  the  lifetime  for 
absorption,  then  the  signature  is  broadened.  The  above  observa¬ 
tional  fact  apparently  shows  that  this  effect  is  also  not 
important  --  a  result  that  is  more  persuasive  than,  but  consistent 

to  an  order  of  magnitude  with,  an  estimate  of  the  time  constant  for 

2 

diffusion,  which  is  r  /4d  «  7  months. 

& 

In  Van  Allen  [1982]  it  is  shown  that  the  energy  distribution 
of  protons  in  the  vicinity  of  Ring  G  is  such  that  about  half  of  the 
particles  have  energy  Ep  greater  than  200  MeV,  and  half  have  energy 


Ep  1-sr  than  200  MeV.  This  is  the  foundation  for  adopting  200  MeV 
as  a  representative  value  of  in  the  foregoing  numerical  examples. 

Overall,  I  estimate  that  the  numerical  values  in  equations 
(lM  and  (l6)  are  trustworthy  to  a  factor  of  two. 
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ADDITIONAL  COMMENTS 

A  noteworthy  feature  of  equation  (13)  is  that,  for  b 
sufficiently  small  (i.e.,  ^  10  tan),  <P>  is  proportional  to  the 
sum  of  the  cross -sections  of  the  assemblage  of  absorbing  satellites 
as  written  in  (15)-  The  latter  equation  is  applicable  to  particulate 
matter  in  the  form  of  small  spheres  of  radius  b  and  volumetric  mass 
density  p  provided  that  4  p  b/3  is  greater  than  the  range  of  the 
relevant  protons  («  26  g  cm  )  in  the  particulate  material.  If 
a  sphere  is  pulverized  into  n  equal  smaller  spheres  and  if  the 
latter  are  dispersed  widely  (no  "shadowing"),  then  the  projected  area 
of  the  assemblage  is  increased  by  the  factor  n1^.  If,  however, 
the  radii  of  the  small  objects  are  such  that  4  p  b/3  is  much  less 
than  the  range  of  the  protons,  a  given  quantity  of  material  becomes 
an  even  more  effective  absorber.  The  latter  case  is  found  to  be 
applicable  to  the  particulate  matter  in  Ring  G,  if  its  absorption 
signature  is  attributed  entirely  to  small  particulates  of  uniform 
size  [Van  Allen,  1983].  The  inferred  sum  of  the  volumes  of  all 
particulates  in  the  ring  is  equivalent  to  the  volume  of  a  single 
sphere  of  radius  only  0.1  tan.  The  inferred  total  number  of 
particulates  is  4.3  x  10^.  By  the  pulverization  theorem,  the 
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radius  of  a  single  satellite  having  the  same  absorption  for 
energetic  protons  would  be  8  1cm  if,  contrary  to  my  previous 
analysis,  the  radii  b  of  the  particulates  were  greater  than 
U  p  b/3 . 

For  a  single  satellite  with  b  >  100  km,  <  >  must  be 

set  equal  to  unity  for  each  of  the  two  opportunities  during  an 
encounter.  Then 


<P>  =  2<p2> 


2rrb 


V  T„  cos  a 
B  o 


(19) 


or  numerically 


<  P  > 


3.5  X  10"5  b 


(20) 


for  aQ  =  70°  and  b  in  km.  In  this  case,  <P>  is  proportional  to 
the  first  power  of  b  and,  corresponding  to  (l6)  and  (17) 


T  ,  19§ 


years 


for  a  single  satellite  or 


T  = 


Z  b. 


years 


(21) 


(22) 


for  mor°  than  one. 


The  intersection  of  (20)  and  (l4)  occurs  at 

b  =  100  km.  ( 

There  is  a  gradual  transition  from  the  quadratic  regime  for 
b  ^  20  km  to  the  linear  regime  for  b  ^  100  km.  As  b  increases 
further,  < P >  approaches  1.0  asymptotically. 


CONCLUSIONS 


The  analysis  in  Van  Allen  [1983]  finds  that  the  lifetime 

against  absorption  of  200  MeV  protons  whose  guiding  centers  lie 

within  ±  r  of  the  center  of  Ring  G  lies  in  the  range  1.6  to 
§ 

3-5  years.  If  one  adopts  an  intermediate  value  of  2-5  years, 
then  by  (l6) 


b  -  3*3  km 


(2k) 


or  by  (17) 


2  2 

I  b  =  11.2  km  .  (25) 

i=l  1 

The  reader  is  reminded  that  these  estimates  are  in  the  nature  of 
extreme  upper  limits  for  the  radii  of  one  or  more  shepherding 
satellites  because  they  ignore  the  absorbing  effects  of  the 
optically-observed  distribution  of  particulate  matter. 

It  is  difficult  to  assign  a  rigorous  uncertainty  to  the 
foregoing  values  but  I  estimate  that  equation  (2k)  is 
uncertain  by  less  than  a  factor  of  two  and  equation  (25'>,  by 
less  than  a  factor  of  four. 
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Table  1.  Pitch-Angle  Dependence  of  the  Probability  of  Absorption 
--  Relative  to  That  at  70°. 
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FIGURE  CAPTIONS 


Figures  1(a)  and  1(b). 

The  origin  of  the  adopted  Cartesian  coordinate  system  is 
the  point  at  which  the  guiding  center  of  a  particular  gyrating 
particle  pierces  the  equatorial  plane;  the  origin  moves  in 
longitude  at  the  same  average  rate  as  does  the  guiding  center. 

The  Y-axis  if  directed  toward  the  planet's  center  and  the  Z-axis 
is  parallel  to  its  axis  of  rotation.  In  each  of  the  two  drawings, 
the  large  circle  represents  the  equatorial  projection  of  the  seg¬ 
ment  of  the  particle's  helical  trajectory  that  lies  near  the 


equatorial  plane.  The  small  circles  represent  successive  positions 
of  the  equatorial  cross-section  of  a  spherical  satellite  as  it 


moves  at  constant  speed  v  parallel  to  the  instantaneous  X-axis 

with  its  center  at  y  =  yQ,  z  =  0.  Note  that  for  the  cases  of 

interest  in  this  paper  b/r  <  0.001.  A  much  larger  value  of 

g 

this  parameter,,  namely  0.1,  was  chosen  for  these  drawings  in 


the  interest  of  graphical  clarity. 


Figure  2 


A  composite  plot  of  equations  (l)  and  (2)  of  the  text 

in  dimensionless  form.  Note  that  this  illustrative  plot  is  for 

b/r  =  0.01,  a  value  of  the  order  of  ten  times  those  of  relevance 
'  g 

to  the  cases  of  interest  in  this  paper. 
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